Orientational Ordering in Sequence-Disordered Liquid Crystalline Polymers 
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Phase separation of sequence-disordered liquid crystalline polymers, a promising class of tech- 
nological and biological relevance, is studied by field theory, and thermodynamic mechanisms re- 
sponsible for orientational ordering observed in experiments, are discussed. The theory developed 
predicts that chemical disorder marginally affects the nematic/isotropic biphasic coexistence width, 
but strongly impacts ordering; above a critical chemical disorder threshold orientational ordering is 
precluded. 
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The study of sequence heterogeneity effects on phase 
separation of sequence-disordered Uquid crystalline poly- 
mers is paramount to modulation of electro-optical prop- 
erties, piezo and pyro-electricity and also to disclo- 
sure of general ordering trends prevalent in proteins with 
secondary structure ||^. Onsager approach devel- 
oped for the study of liquid crystals (LC's) with repulsive 
anisotropic interactions, has been extended and used ^] 
to analyze single chains and melts of liquid crystalline ho- 
mopolymers (TCP's). The effect of attractive anisotropic 
interactions by dipolar coupling was first introduced in 
the theory of TC's by Mayer and Saupe ||^; their ap- 
proach was found useful in the study of low molecular 
weight TC's. TCP's were also studied by lattice models 
and the Onsager approach was used to correct for 
lattice artifacts effects in phase behavior 0. 

A field theoretic approach to study many chain TCP's 
was proposed by Gupta and Edwards (GE) free 
energy contributions to nematic ordering were com- 
puted non-perturbatively, and the Tandau expansion 
method was avoided by summation over all Rouse modes; 
interaction-wise, both athermal repulsive contributions 
1^, and the Mayer Saupe dipole-dipole attractive con- 
tributions to the inter-segment interaction potentials are 
present. GE |H predict that the nematic/isotropic (N/I) 
transition is first order while the transition temperature 
is depressed in comparison with low molecular weight 
rods; their predictions are supported by computer simu- 
lations of short nematic chains pT[ | and by experiments 
. Recently, GE approach was employed to study ori- 
entational ordering in semi-flexible homopolymers em- 
bedded in flexible surfaces |12|. 

Stupp et al. 1^, synthesized and characterized 
by optical microscopy and C13 NMR, orientation and 
sequence-statistics of thermo-tropic sequence-disordered 
TCP. In these experiments, optical domains were ob- 
served over a broad temperature regime '-^ 120°; these do- 
mains were attributed to a wide biphasic N/I coexistence 
width, and their occurrence was attributed to sequence- 
disorder and finite chain length effects. Fredrickson et al. 
1^^ , studied via Tandau theory the homogeneous chain 
anisotropy limit of the sequence model , and predicted 
that an increase in chain length diminishes chemical het- 



erogeneity effects on the N/I biphasic width. 

In order to make closer connection with experiments, 
in the present work we construct a field theory of TCPs 
made of stiff mesogens and flexible spacers, randomly dis- 
tributed on the chain (viz. fig. 1). Sequence-distribution 
and interaction-wise, our model is an adequate descrip- 
tion of the sequence disordered TCPs synthesized in ex- 
periments ||T^. While the GE method is not applicable 
here, in the homopolymer limit the present study repro- 
duces the free energy of GE theory 

The many-chain Hamiltonian for a solution of se- 
quence disordered TCP's made of mesogenic segments 
(A's) and flexible segments (B's) is: 
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r(ni) is the spatial location of the n'th segment of the 
i'th chain, u(ni) is the chain tangent of the i'th chain 
at ni and 0{ni) is the chemical composition variable 
of the n^'th segment on the i'th chain. 9{ni) = 1 
for an A segment and 0{ni) = —1 for a B segment. 
The sequence heterogeneity, represented by fluctuations 
in the segment composition along the chain contour 
obeys a Gaussian process with mean, < 9 >— 2/ — 1, 
S9 = 9{ni)— < 9 >, and sequence fluctuations 69'^ =< 
iS9{n,)69{n'^)) >=6{n, - n;-)4/(l - /)/ [|l§; f is the frac- 
tion of A segments, and 1 is the statistical segment length. 
The first term in eq. ^ is the nearest- neighbors harmonic 
interaction potential of flexible segments of type B, while 
the second term in eq. Q precludes bending of the stiff 
segments, A's, (viz. line-bounded mesogens, fig. 1). The 
third term in eq. |^ is the interaction potential of A- A seg- 



this energetic penalty is zero for aligned tangents of A 
segments adjacent in space, and w for normal tangents. 




Fig. 1: One disordered LCP sequence. Flexible segments 
are thin and black, mesogen segments are thick and grey 

The A-A anisotropic interaction potential, w, con- 
tains both athcrmal and thcrmo-tropic contributions . 
The fourth term in eq. 0, is the excluded volume inter- 
segment interaction. While sequence effects on equi-stiff 
heteropolymers were treated the present work fo- 
cuses on the effect of dissimilarity in anisotropy among 
the stiff mesogens (A) and the flexible spacer (B) on 
orientational ordering, in lyotropic and thermo-tropic 
sequence-disordered LCPs. 

Below, solution steps are briefly described. First, con- 
tinuous microscopic orientational tensors are introduced 
by delta function constraints on the partition function, 
while the fluctuations of the chain director of the stiff 
segments is constrained in the usual way (viz. [|j for fur- 
ther details). The partition function at fixed sequence is 
given by: 

Z[0K),u(nO]= J Jill Dr,{n,)D<7{r)]exp{-L[a{r)]) 



tational convenience, we exchange the Lagrangian to 
a Hamiltonian representation and the free energy 
is computed variationally in the Hamiltonian reference 
frame: 



Ha = '^{-JJPI + haul) 
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and are the a quantum kinetic energy and coor- 
dinate operator in imaginary time with a = x,y, z. 4'a 
is the principal axis representation of the fields y/—lip'^^ 
while A equals ^/—l times the auxiliary field that sets 
the magnitude of u{ni) to 1 for the A segments in eq. ||. 
Based on eq. ^, creation/annihilation operators, (a"*", a), 
are introduced to facilitate computation of orientational 
averages in the free energy. 
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Omitting presentation of the intermediate calculation 
p2|, the free energy per segment obtained is: 
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p is the segment density; a a is (Ji.j in principal axis rep- 
resentation: 



a{r) is given in eq. ^; the d'(r) delta function constraint 
is expressed with auxiliary fields ip"^^ (r) , and the disorder 
average is performed non-perturbatively with the replica 
trick pO|] ; replica symmetry breaking (RSB) is not ex- 
pected and the replica limit is exact; still the disorder 
average yields a free energy functional non-local in poly- 
mer modes, and the GE approach [|| cannot be used. 

The partition function in eq. ^ may be also viewed 
as the path integral for a quantum Green function in 



a-b 
cr" = I a + b 
2a 



(8) 



For uniaxial ordering the orientational order parameter, 

< 5 >, is given by < >= 1 > < 5* > > 
signals uniaxial nematic ordering while —0.5 < < S > 

< signals discotic ordering. In Einstein notation, in 
the isotropic phase, a"{r) = pa(i"), otherwise, fT**(r) 7^ 
PA(r). The homopolymer limit of the free energy, i.e. 



a many-chains LCPs (viz. eq. 32 in Q), and the min- 
imal value for uniaxial ordering, < S >=0.25, obtained 
analytically in S), is reproduced herein. 
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Fig. 2: Numerical results for uniaxial ordering < 5 > 
(solid lines) and chemical potential per segment (long dashed 
lines) variation with segment density, p. Comparison among 
the thermodynamical stable phase (thick line) and the unsta- 
ble phase (thin line); f=0.938, w=10, /3e=5, 1=8. 

Now, the free energy given in eq. |^ is minimized 
analytically by the fields aa,ipa,^ with a — x,y,z] 
seven non-linear self-consistent equations are obtained 
and solved numerically by direct iteration with proper 
mixing of coefficients ^] ; more computational and ana- 
lytical details will be given elsewhere |^ . Fig. 2. depicts 
numerical results for lyotropic stability of disordered LCP 
at small disorders. Thick lines in fig. 2 are uniaxial order- 
ing curves for < S" > while dotted lines are the chemical 
potentials of the stable phases. 




Fig. 3: Numerical results for uniaxial ordering < 5 > 
variation with segment density, p and composition, f; w=10, 
/3e=5, 1=8; f values follow upper to lower < S > curves; 
f=0.98, 0.94, 0.93, 0.88. 

T'Vici f.Viir*k linci aVir»isriTi(r q rliGr*r»nf iTinr^na inmn in ,^ 



at a finite density depicts a first order N/I transition; fig. 

2 depicts also another stable solution for < S" > (viz. 
other thick line); this solution increases continuously with 
increase in density from < S >=0, i.e. a second or higher 
order N/I transition. 

The first order line bears a lower chemical potential 
at all densities studied (viz. thick dashed line), and de- 
scribes the physical ordering scenario that should be ob- 
served in simulation or experiments. The weak N/I tran- 
sition, is also numerically stable, but has a higher chemi- 
cal potential (viz. thin dashed line) and it should not be 
observed. 

Fig 3. is a numerical study of sequence heterogeneity 
effects on lyotropic nematic ordering. Only stable solu- 
tions with lowest chemical potentials are plotted. Fig. 

3 depicts < 5 > - the nematic order parameter varia- 
tion with density for decreasing fractions of the stiff seg- 
ments, f. The entropy carried by the flexible runs have 
dramatic effects on ordering; Fig. 3 shows that a rel- 
atively small increase in the average number of flexible 
segments decreases notably the overall uniaxial ordering, 
and shifts the segment density at the N/I threshold to 
higher values. Interestingly, while the presence of se- 
quence heterogeneity impacts strongly the value of the 
segment-density-N/I-threshold, the presence of disorder 
carried by the fluid does not change the order of the N/I 
transition, and the transition remains flrst order. 
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Fig. 4: Numerical results for N/I biphasic width variation 
with fraction of mesogen segments, f; w=10, /3e=5, 1=8. 

Recently Radzihovsky et. al. have shown that 
weak external disorder is sufficient to destabilize the ne- 
matic phase in LC rods; in the present system, wherein 
the disorder is carried by the fluid, nematic ordering is 
also affected; our calculation show that above a critical 
value of disorder strength, the nematic ordering is de- 
stroyed and the phase becomes isotropic, yet the overall 
effect on uniaxial ordering is less pronounced then the 
effect of external disorder. 

A Maxwell construction is employed (viz. [0]) for the 



calculation of the N/I coexistence region. Free energies 
and chemical potentials in the vicinity of the N/I tran- 
sition are computed numerically for the isotropic and 
the nematic phases for several fractions of mesogenic 
segments. A, the N/I coexistence width is given by 
A = Pi — pn] Pi is the segment density in the isotropic 
phase while pn is the segment density in the nematic 
phase at coexistence. Fig. 4 is a study of sequence het- 
erogeneity effects on the N/I coexistence. In fig. 4 nu- 
merical results for the coexistence width dependence on 
mesogenic fraction are depicted. Larger sequence het- 
erogeneities increase the N/I coexistence width almost 3 
folds that of a semi- flexible homopolymer. 
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Fig. 5; N/I phase diagram in f and p variables; w=10, 
/3e=5, 1=8; squares are {i,pi) points, circles are (f,pn) points. 

Fig. 5 depicts the effect of sequence heterogene- 
ity on orientational ordering. The entropy carried by 
mesophases rich in flexible segments destroys the nematic 
ordering below a specific fraction of mesogen segments 
(which is density dependent), the nematic phase loses sta- 
bility, while the isotropic phase remains stable. For the 
values of parameters in fig. 5 at melt densities, {p — 1), 
the critical fraction of mesogen segments, f, is fc— 0.76. 
The N/I coexistence width displayed in fig. 5 as the hor- 
izontal distance among an adjacent square (pi) and circle 
(pn) is marginal at all densities displayed. 

Invoking the homogeneous isotropic limit, (i.e. the 
difference in anisotropy of A and B segments is negligible) 
the sequence heterogeneity effects vanish in the infinite 
chain limit in agreement with |p^ ; In many experimen- 
tal scenarios, (i.e. [MBPE — {methylene)y with y=17, 
18, 20 ll^) one segment is stiff, the other is flexible and 
rigorous treatment of anisotropy dissimilarity of A and B 
segments developed herein is crucial. In the present work, 
we showed that for long chains, ^/L << 1, the sequence 
heterogeneity effect on the N/I biphasic region does not 
vanish, but in most cases is small with a AT ~ 10°. 



contribution to the N/I biphasic width. Our predictions 
suggest that optical mesophases in the regime ^/L << 1 
cannot be associated with a wide N/I biphasic regime in- 
duced by sequence-heterogeneity, and a different expla- 
nation of experiments, based on a rigorous account of the 
optical scale ordering of mesophases (e.g. cr(k) ^ 0), is 
necessary. The occurrence of optical domains observed in 
experiments, |l5| , can also be a signature of slow order- 
ing dynamics during a finite heating/cooling rate process. 
Currently, we are studying both scenarios. 
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